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We develop a theory of nonlinear cosmological perturbations on superhorizon scales for a multi- 
component scalar field with a general kinetic term and a general form of the potential in the context 
of inflationary cosmology. We employ the ADM formalism and the spatial gradient expansion 
approach, characterised by 0{e^), where e = 1/{HL) is a small parameter representing the ratio 
of the Hubble radius to the characteristic length scale L of perturbations. We provide a formalism 
to obtain the solution in the multi-field case. This formalism can be applied to the superhorizon 
evolution of a primordial non-Gaussianity beyond the so-called 5N formalism which is equivalent to 
©{e") of the gradient expansion. In doing so, we also derive fully nonlinear gauge transformation 
rules valid through C(e^). These fully nonlinear gauge transformation rules can be used to derive 
the solution in a desired gauge from the one in a gauge where computations are much simpler. As 
a demonstration, we consider an analytically solvable model and construct the solution explicitly. 

PACS numbers: 98.80.-k, 98.90.Cq 



I. INTRODUCTION 

■ Recent observations of the cosmic microwave background anisotropy show a very good agreement of the observational 
'"^ data with the predictions of conventional, single-field slow-roll models of inflation, that is, adiabatic Gaussian random 
primordial fluctuations with an almost scale-invariant spectrum Nevertheless, possible non-Gaussianities from 
inflation has been a focus of much attention in recent years, mainly driven by recent advances in cosmological 
observations. In particular, the PLANCK satellite 0] is expected to bring us preciser data and it is hoped that 
a small but finite primordial non-Gaussianity may actually be detected. 

To study possible origins of non-Gaussianity, one must go beyond the linear perturbation theory. An observationally 
detectable level of non-Gaussianity cannot be produced in the conventional, single-field slow-roll models of inflation, 
since the predicted magnitude is extremely small, suppressed by the slow-roll parameters. Then a variety of ways to 
generate a large non-Gaussianity have been proposed. (See e.g. a focus section in CQG Q and references therein 
for recent developments.) They may be roughly classified into two; multi-field models where non-Gaussianity can 
be produced classically on superhorizon scales, and non-canonical kinetic term models where non-Gaussianity can be 
produced quantum mechanically on subhorizon scales. In particular, in the former case, the 6N formalism turned out 
to be a powerful tool for computing non-Gaussianities thanks to its full non-linear nature. 

On the superhorizon scales, one can employ the spatial gradient expansion approach P-[30j. It is characterised 
by an expansion parameter, e = 1/{HL), representing the ratio of the Hubble radius to the characteristic length 
scale L of the perturbation. In the context of inflation, based on the leading order in gradient expansion, the 6N 
formalism [1, [ill, [Hj or the separate universe approach [l^ was developed. It is valid when local values of the infiaton 
field at each local point (averaged over each horizon-size region) determine the evolution of the universe at each point. 
This leading order in the gradient expansion provides a general conclusion for the evolution on superhorizon scales 
that the adiabatic growing mode is conserved on the comoving hypcrsurfacc [l7| . 

In this paper, we consider the curvature perturbation on superhorizon scales up through ncxt-to-lcading order in 
gradient expansion, that is, to ©(e^). To make our analysis as general as possible, we extend the SN formalism in the 
following two aspects: One is to go beyond the single-field assumption, and the other is to go beyond the slow-roll 
condition. While in the case of single-field inflation, the curvature perturbation remains constant as mentioned above, 
the superhorizon curvature perturbation can change in time in the case of multi-field inflation. Furthermore, even for 
single-field inflation, the time evolution can be non-negligible due to a temporal violation of the slow-roll condition. 
In order to study such a case, the SN formalism is not sufficient since the decaying mode cannot be neglected any 
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longer, which usually appears at Qie^) of gradient expansion and is known to play a crucial role already at the level 
of linear perturbation theory [Sll, [32|/ not to mention the case of nonlinear perturbation theory [l^, HH, [28l. [29|. 

Multi-field inflation may be motivated in the context of supergravity since it suggests the existence of many flat 
directions in the scalar field potential. In multi-field inflation, a non-slow-roll stage may appear when there is a change 
in the dominating component of the scalar field. For example, one can consider a double inflation model in which a 
heavier component dominates the first stage of inflation but damps out when the Hubble parameter becomes smaller 
than the mass, while a lighter component is negligible at the first stage but dominates the second stage of inflation 
after the heavier component has decayed out [24. 133. [ssj . However, these previous analyses are essentially based on 
the 5N formalism and it is in general necessary to extend it to 0{e'^), that is, to the beyond 6N formalism. We focus 
on the case of a multi-component scalar field. As for a single scalar field, it has been developed in 12M . 

We mention that a multi-scalar case in the gradient expansion approach was studied previously |25l. [26j. However, 
it turns out to be valid only for a restricted situation (discussed later) . Here we develop a general framework for fully 
nonlinear perturbations and present a formalism for obtaining the solution to 0{e^). Then as an example we consider 
a specific model which allows an analytical treatment of the equations of motion. 

This paper is organised as follows. In Sec. [Hi we introduce a multi-component scalar field and derive basic equations. 
We compare several typical time-slicing conditions and mention the differences of them from the single-field case. In 
Sec. mil we develop a theory of nonlinear cosmological perturbations on superhorizon scales. We formulate it on 
the uniform e-folding slicing (which is defined later). In Sec. IIVI as a demonstration of our formalism, we consider 
an analytically solvable model and give the solution explicitly. Sec. |V] is devoted to a summary and discussions. 
Some details are deferred to Appendices. In Appendix \^ the coincidence between some of time slicing conditions 
is discussed by using the Einstein equations. In Appendix |Bl we write down the basic equations on the uniform 
expansion slicing, and study the behaviour of the curvature perturbation in this slicing. In Appendix [Cl we give 
general nonlinear gauge transformation rules valid to next-to- leading order in gradient expansion. In Appendix |D1 wc 
verify our formalism in a single-field model. Finally, in Appendix [El we discuss the structure of the Hamiltonian and 
momentum constraint equations in the gradient expansion. 



II. BASICS 



A. The Einstein equations 



We develop a theory of nonlinear cosmological perturbations on superhorizon scales. For this purpose we employ 
the ADM formalism and the gradient expansion approach. In the ADM decomposition, the metric is expressed as 

ds^ = gf,^dx^'dx'' = -a^dt^ + % {dx' + l3'dt) {dx^ + (3^ dt) , (2.1) 

where a is the lapse function, is the shift vector and Latin indices run over 1,2 and 3. We introduce the extrinsic 
curvature Kij defined by 

K^, = ^ {da,j - A/3, - A a) . (2-2) 

where D is the covariant derivative with respect to the spatial metric 7^ j . In addition to the standard ADM decom- 
position, the spatial metric and the extrinsic curvature are further decomposed so as to separate trace and trace-free 
parts 

7,, = (t)e''^7.j ; dct 7^ = I , (2.3) 

^ a\t)e''^ (l-Kj,, + Aj) ; - , (2.4) 



where a{t) is the scale factor of a fiducial homogeneous Friedmann-Lemaitre- Robertson- Walker (FLRW) spacetime 
and the determinant of jij is normalised to be unity and Aij is trace free. The explicit form of K is given by 

K = f^K,, = i [3{H + dt^) - A/3*] , (2.5) 



See, however, a special case of single-field infiation studied recently in [sSl where the would-bc decaying mode of the comoving curvature 
perturbation happens to be rapidly growing outside the horizon and an extended version of the 5N formalism remains to be valid 
although the curvature perturbation is no longer conserved. 
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As for a matter field, let us focus on a minimally-coupled multi-component scalar field, 



where H is the Hubble parameter defined by H{t) 



(2.6) 



where /, J and K run over 1, 2, • • • , with 0^ denoting the K-\h component of the scalar field. Note that we do 
not assume a specific form of both the kinetic term and potential, which are arbitrary functions of X^'^ and 0^. This 
type of Lagrangian can be applied to, for example, multi-field DBI inflationary models. For the calculation of their 
non-Gaussianities, see, e.g. [Sy, H^l and also [H, for recent developments. 
The equation of motion for the scalar field is given by 



where the subscript / in Pj represents a derivative with respect to (f>^ and P{ij) is defined as 



P 



dP 



dP 



[1.1) 



2 \dXiJ dXJi 



The energy-momentum tensor is 



(2.7) 



(2.8) 



(2.9) 



Notice that this energy-momentum tensor cannot be written in the perfect fluid form any more, which is one of main 
differences from the single-field case. 

All the independent components of the energy-momentum tensor are conveniently expressed in terms of E and Ji 

as 



and Tij , where n'* is the unit vector normal to the time constant surfaces and is given by 



(2.10) 



(2.11) 



For convenience, we further decompose in the same way as Eq. ()2.4p . 



a^(t)e 



2-4) 



1 



(2.12) 



Now we write down the Einstein equations. In the ADM decomposition, the Einstein equations are separated into 
four constraints, the Hamiltonian constraint and three momentum constraints, and six dynamical equations for the 
spatial metric. The constraints are 



R - [AD^i) + 2DVA^/'jJ + gA'- 



'-A,jA^^ = 2E, 
(^^^ A-^ — Ji , 



(2.13) 
(2.14) 



where R = R['y\ is the Ricci scalar of the normalised spatial metric 7^-/, Di is the covariant derivative with respect to 
7jj, = j'-^DiDj, 7*-' is the inverse of 7ij, and the spatial indices are raised or lowered by 7*-' and jij, respectively. 
As for the dynamical equations for the spatial metric, we rewrite Eq. (|2.2p as 



d±il> = + - A" + 

a 3 V ct 



7»j 



= 2A,, 



TF 



(2.15) 
(2.16) 
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The equations for the extrinsic curvature (K^Aij) are given by 

d±K = -1^2 _ ^ + ^_fD^a + D'aD.i) -\{S + E) , (2.17) 



1 



2p2'i/j 



1 ' -'jn.-'ir' Q 

ay o 

i?y + Dii^Dji! - D^Djtjj - - [ DiDja - DiaDji/j - DjipDia 



1 TF 



S^J , (2.18) 



where d± — n^d^, and we have introduced the trace-free projection operator [...J"'"^ defined for a tensor Qij as 



QV = Qv - llvl'^Qki ■ (2.19) 



Finally, the equations of motion for the scalar field (|2.7p are 



d 



{P(ij)d^(t>')+KP(ij)d^<l,' - -l_.d,(aae^P^ij)f^dj<ly') -\Pi^^- (2-20) 



B. Gradient expansion and assumption 

In the gradient expansion approach we suppose that the characteristic length scale L of a perturbation is longer 
than the Hubble length scale \/H of the background, i.e. HL ^ 1. Therefore, e = 1/{HL) is regarded as a small 
parameter and we can systematically expand equations in the order of e, identifying a spatial derivative is of order e, 
diQ = 0{e)Q. To clarify the order of gradient expansion, we introduce the superscript (n). For example, *^^)q; means 
the lapse function at second order in gradient expansion. 

As a background spacetime, we consider a FLRW universe. At 0(6°) of the gradient expansion, there is apparently 
no spatial gradient and the universe is locally homogeneous and isotropic. This leads to the following condition on 
the spatial metric: 

^a^J - O(e^) . (2.21) 

Since we adopt this assumption, the spatial metric at leading order is given by an arbitrary spatial function of the 
spatial coordinates, 

'•'H.^hi^'), (2.22) 

under the condition that the eigenvalues of fij are all positive definite everywhere. From the definition of Aij, 
Eq. (|2.16p . the above assumption implies 

Ay = 0{e^) . (2.23) 

Throughout this paper, in order to simplify the equations, we set the shift vector to zero up to second order in 
gradient expansion, 

= Oie^) . (2.24) 

Let us call this choice of the spatial coordinates as the time- slice- orthogonal threading. Here we mention that the 
above condition does not completely fix the spatial coordinates. As discussed later, one can actually make an arbitrary 
coordinate transformation of the form, — > = /'(x^ ). 

C. Leading order in gradient expansion 



In this subsection, we study the leading order gradient expansion and make clear the correspondence between the 
leading order equations and background equations. This correspondence can be used to construct the solution at 
leading order in gradient expansion in terms of the background solution. 
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At leading order in gradient expansion, the Einstein equations are 



d,K = -3P(,j)a,0^9,0'^ , 

drK=-iPnj)dr¥drCt>\ 



and the scalar field equation is 



dr (P(IJ)dA') + KP^jj^dr<f>'' - = , 

where we have introduced the proper time r by 



ait', x') dt' . 



(2.25) 

(2.26) 
(2.27) 

(2.28) 
(2.29) 



—const. 



In terms of r, the expression of K in Eq. (12. 5p is simplified under the time-slice-orthogonal threading condition, 



Under the identifications. 



K = \ ' = 3 \ ^ 



ae^ a , and t t , 



(2.30) 



(2.31) 



one also has the correspondence, K 3H. This means that the basic equations at leading order, Eqs. (j2.27p 
and (|2.28p . take exactly the same form as those in the background modulo above identifications. Namely, given a 
background solution. 



^BG 



background 

one can construct the solution at leading order in gradient expansion as 



gradient 



(2.32) 



(2.33) 



All the information of inhomogeneities is contained in the initial condition as well as in the proper time r through 
Eq. (|2.29p . Thus it is sufficient to solve the background equations to obtain the solution at leading order in gradient 
expansion. 

In passing, we note that the e-folding number is often used as the time coordinate to describe the background 
evolution. For convenience, we define it as the number of e-folds counted backward in time from a fixed final time. 
That is. 



N{t) 



to 



H{t') dt' . 



(2.34) 



Accordingly, the scale factor is expressed as 

a{N) = aoe-^^-^") . (2.35) 
By replacing t with r and H with A'/3 we can generalise the e-fold number to the one defined locally in space as 



N{t,x') 



3 



to 



dt' a{t',x')K{t',x') 



(2.36) 



Again one can check the validity of the above correspondence by rewriting Eqs. (|2.25p . (|2.27p and (|2.28p in terms of 
Af as the time coordinate. 
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D. Various slicings and their coincidences 

One needs to specify the gauge condition to study perturbations in perturbation theory or in gradient expansion. 
Since spatial coordinates have been aheady fixed by the time-shce-orthogonal threading, one has to determine the 
timc-shcing condition. Here, let us list various slicings and their definitions, 

Comoving ; J^ = , (2.37) 

Uniform expansion ; K{t,x') = 3H{t) , (2.38) 

Uniform energy ; E{t, x') = E{t) , (2.39) 

Synchronous ; alt.x^) = I , (2.40) 

Uniform e-folding number ; 7V(t, a;*) = N{t) . (2-41) 

Hereinafter, we call the uniform expansion, uniform energy and uniform e-folding number slicings as the uniform K, 
uniform E and uniform Af slicings, respectively. 

We mention that there is a remaining gauge degree of freedom in the synchronous or uniform Af slicing, while the 
time slices are completely fixed in the uniform K and uniform E slicings. As for the uniform J\f slicing, the gauge 
condition demands dtip to vanish from Eqs. (|2.30p and (|2.36p . This means one can freely choose the initial value of 
"0 (and hence its spatial configuration at any later time because "0 is conserved). This corresponds to the freedom in 
the choice of the initial time-slice as we see later. Utilising this freedom, we can make a scalar quantity, one of scalar 
fields 0^ or K for example, homogeneous on the initial slice. 



E. Towards the next-to-leading order in gradient expansion 

As we have seen in subsection lH CI the leading order solutions are given by functions of r in terms of the background 
solutions. At ncxt-to-lcading order in gradient expansion, terms with spatial derivatives of the leading order solution 
appear in the evolution equations. To evaluate those terms, one needs to calculate the spatial derivative of the lapse 
function, for example in dicf), 



d.,^BG{T) = a,0BG(r) a,T = a,0BG(r) J d^adt . (2.42) 
However the leading order ^^"^a is in general given explicitly only after solving the following equation for a: 



f t, (j>{T) 



f 



t, 4> [ adt 



(2.43) 



As a demonstration, the analysis on the uniform K slices is performed in Appendix [Bl and it is clearly shown that it 
is almost impossible to solve this equation, at least in an analytical way. 

This problem did not appear in the single-field case. It is because one can show that various different slicings become 
identical at leading order in gradient expansion. In particular, all the slicings listed in subsection III Dl coincide with 
each other as shown in Appendix [X] 

comoving = uniform K = uniform E <--- synchronous = uniform J\f , (2.44) 

where — >■ means the left ones imply the right ones and <--- means it holds when one chooses the initial slice to be the 
comoving, uniform K or uniform E slicing by using the remaining gauge degree of freedom. Thus the lapse function 
is homogeneous in all the slicings in the above, and we may set a = I if desired. 

On the other hand, one has to face this problem in the case of multi-field inflation. We overcome this problem 
by choosing the synchronous slicing or uniform M slicing, which gives us a homogeneous time coordinate. On these 
slicings, one can evaluate the spatial derivatives of the leading order solution which appear as source terms and 
construct a solution to next-to-leading order in gradient expansion by integrating those terms. 

There are two necessary steps before reaching the goal. Once we have a solution, it is necessary to construct a 
conserved quantity out of it that can be directly related to observable quantities. It is widely known that the comoving 
curvature perturbation eventually become conserved in a single-field model in linear theory. In non-linear theory, there 
exists a corresponding quantity, ip on the comoving, uniform K or uniform E slicing, which is conserved at leading 
order in gradient expansion [17[. Even in the multi- field case, the system effectively reduces to a single-field system 
after the so-called non-adiabatic pressure has died out, that is, when the adiabatic limit is reached. Therefore it is 
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necessary to perform a nonlinear gauge transformation from the uniform slicing to one of those three slicings. This 
is one of the steps. Since the comoving slicing is not well defined in general in the multi- field case [l2|, we choose the 
uniform K slicing as the target gauge. 

The other step to be taken is related to the definition of the curvature perturbation at next-to-leading order. In 
linear theory the curvature perturbation is named so because it determines the three-dimensional Ricci scalar, and ip 
can be called so to full nonlinear order in the context of the leading order in gradient expansion. At next-to-leading 
order, however. ■0 itself is no longer adequate to be called the curvature perturbation [28| . One needs to add the 
contribution from part of 7^, which we call Xj to obtain a properly defined curvature perturbation conserved through 
0{e^). Therefore, after transforming from the uniform M slicing to the uniform K slicing, one has to evaluate the 
combination, ^Kk = V'if + Xk/"^- This is the Beyond SN formalism. 

Before concluding this section, we mention the difference between our work and that of Weinberg [2^, [2^ . There it 
was assumed that the lapse function can be chosen to be equal to unity at leading order in gradient expansion, hence 
all the scalar fields are homogeneous. This severely constrains the class of scalar field models as well as the initial 
condition because the curvature perturbation must be always conserved at leading order in gradient expansion. Here 
we do not impose such assumptions and perform a completely general analysis. 



Let us first summarise the five steps in the Beyond SN formalism. 

1. Write down the basic equations (the Einstein equations and scalar field equation) in the uniform M slicing 
with the time-slicc-orthogonal threading. For convenience let us call the choice of the coordinates in which one 
adopts the uniform X slicing with the time-slice-orthogonal threading the X gauge. So the above choice is the 
Af gauge. In this gauge the metric components at leading order are trivial since both ip and 7,,, are independent 
of time. 

2. First solve the leading order scalar field equation under an appropriate initial condition and then the next-to- 
leading order scalar field equation which involves spatial gradients of the leading order solution. 

3. Solve the next-to-leading order Einstein equations for the metric components and their derivatives. 

4. Determine the gauge transformation from the Af gauge to the K gauge and apply the gauge transformation 
rules to obtain the solution in the K gauge. 

5. Evaluate the curvature perturbation 9\ ^ tp + x/3 in the K gauge, where x is to be extracted from 7^ . 

In what follows, we describe these steps in detail but only formally. An example in which these steps can be computed 
analytically will be discussed in Sec. IIVI 



III. BEYOND SN FORMALISM 



Step 1: 



First, we rewrite the uniform slicing condition from Eqs. (|2.30l) and (I2.36P as 



a{t,x')K{t,x') =3H{t) ^ dt^p{t,x') =0. 



(3.1) 



Hence ip is constant in time and is given by a function of the spatial coordinates alone. 



^{t,x') = ^ix') EE C"^(a;*). 



(3.2) 



In the Af gauge, the Einstein equations are reduced to the following equations. The constraints are 



1 




(3.3) 



9 




(3.4) 
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The evolution equations for K, Aij and are 

3 



1 



1 



dNAij = 2,Ai 



TF 



3 



The scalar field equation is 



/v 



K 



P{ij)dN(, 



P{ij)dN(, 



K 



„c* 



We rewrite Eq. by eliminating di^K with Eq. p.Sp as 



Pnj)l''d,^-'\ --Pi^O. 



(3.5) 

(3.6) 
(3.7) 

(3.8) 



9 



if 



AC"'' 



1 



P{ij)dN(t>'' ■ 



(3.9) 



Thus once is expressed in terms of the scalar field and its derivatives, the above equation gives a closed scalar 
field equation. An explicit derivation of the closed scalar field equation is possible only after we specify the explicit 
form of P{X^'^ , 0) as a function of X^^^ and cj). Here we describe generally but formally the procedure to obtain 
as a function of the scalar field. 

First, we separate the term with time-derivatives in X^"^ and denote it by K^Y^^ where Y^'^ = dN4>^ dN4>^ 



1 



(3.10) 



At leading order, we can neglect the spatial-derivative term in Eq. p.lOp . Then, P and P{ij) are given by 
P{K'^Y^'^ and P(^jj'f{K'^Y^'^ ,(t)^), respectively. To next-to-leading order, P and P(ij) can be expanded as 



P{X",<t>^) = P{K'Y",<p^) - -^Y^d.^'^U'd/'^U''^ +0{e') , 



(3.11) 
(3.12) 



Inserting these expressions into Eq. p.3|) . one obtains an algebraic equation for K^. Solving it gives an expression of 
in terms of the scalar field. Then a closed equation for the scalar field is obtained by plugging it into Eq. p.9p . 

Step 2: 

Although we can keep our discussion completely general, below we focus on the case of a multi-component canonical 
scalar field, 



P=hjjX'-' -V{<j>\--- ,^^) and P^jj^ = lsjj. 



(3.13) 



This choice is taken purely for the sake of simplicity and clarity, because the expansion K can be explicitly expressed 
in terms of the scalar field in this case. In general, one cannot obtain an explicit expression of K in terms of the 
scalar field unless the form of P is explicitly specified. Nevertheless, the discussion below also applies to the general 
case perfectly. 
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From Eq. (jOl) wc find 
/l2 / 1 



9 



1 



6 



dN(t>idN(, 



V 1 

"3 " 6a2e2C'* 



(3.14) 



Inserting this into Eq. p.9p . one obtains the following closed equation: 




3 / e<^* 



V 



dN(j)i } - 



R - (aD^C''' + 2D'C'f'DiC'^^ + D'(j>jD, 



-1 



dl 



2d 



'N9l 



(3.15) 



Since each term in the right-hand side of the equation involves two spatial derivatives, we can neglect them at leading 
order. At next-to-leading order, they can be understood as source terms whose time evolution have been already 
determined from the leading order solution. As noted in the above, although one cannot obtain an equation like the 
above explicitly for general P, one can always derive a closed equation for the scalar field once a specific form of P is 
given. 

Solving the closed equation for the scalar field, the solution is formally obtained as 



(3.16) 



where the subscript Af indicates the Af gauge. In the arguments on the right-hand side, the subscript denotes the 
initial value, and D{- ■ ■ ) the spatial derivatives of the quantities inside the parentheses. 

Step 3: 

Once the solution for the scalar field is obtained, one can solve Eqs. (|3.5|) . ()3.6p and p.7p to obtain the metric 
quantities. As for K, however, it is simpler and in fact better to use Eq. (|3.14p . which is essentially the Hamiltonian 
constraint, since the integral constants appearing from integrating Eq. (|3.5|) are not freely specifiable but must satisfy 
the Hamiltonian constraint. 

Step 4: 

Given the solution in the Af gauge, the next step is to find the gauge transformation from the Af gauge to the K 
gauge. It can be determined as follows. As noted above, the expression for Kj^ in terms of the scalar field is obtained 
from Eq. p.l4p . Since the leading order and next-to-leading order scalar field solutions are expressed as Eqs. (|3.16p . 
the same is true for K^f, 



N; (2)00, ^'^dNcbo, Di^'Uo, (°'a^0o, ("'t^.o) 



(3.17) 



Let the transformation from the uniform Af slicing to the uniform K slicing be given by N ^ N = N + n{N, x*) or 
conversely N + n{N, i') = A^, where the uniform K slice is given by K =const.. This nonlinear gauge transformation 
is discussed in detail in Appendix [Cl The nonlinear gauge transformation generator n{N,x^) from the Af gauge to K 
gauge is then determined by the condition that K is spatially homogeneous on the uniform K slice by definition: 



^°^KKiN, i') = (°)ifA-(A + x') = , 



0. 



(3.18) 
(3.19) 



Then the other quantities in the K gauge are obtained by applying the above gauge transformation. In particular, 
the determinant of the spatial metric tpK in the K gauge is obtained from Eqs. (|C28p and (|C33p . and the unimodular 
part of the spatial metric 'jij k from Eqs. (jC29p and (jC34p . 
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step 5: 

We are to construct the nonlinear curvature perturbation, D\k = V'-ftr + Xif/S, where is the scalar part of the 
metric ^ij k- The scalar part of ^ij is defined in the same way as the single-scalar case |27| . 



(3.20) 



where A ^ is the inverse Laplacian operator on the flat background. Extracting xk from jijK, and combining ipK 
and Xk, wc finally obtain the nonlinear curvature perturbation ~ ipK + Xif/S m the K gauge. 

IV. SOLVABLE EXAMPLE 

In this section, we demonstrate how to obtain the solution up to next-to-leading order in gradient expansion by 
applying our formalism to a specific, analytically solvable model. 



A. Model and equations 



For simplicity, we consider a canonical scalar field with exponential potential (23|, 



where is a constant. The leading order scalar field equation is given by 



(4.1) 



(4.2) 



where we have omitted a summation symbol over the field component indices, J. Hereafter summation is implied 
over repeated component indices. 

Further we assume the two slow-roll conditions on the leading order equation. The first one is that wc can neglect 
the "kinetic energy" in the energy density of the scalar field. 



The second one is that we can neglect the "acceleration" , 



(4.3) 



(4.4) 



It is important that we apply these assumptions only to ''°'(/). Wc do not impose the slow- roll conditions on (2)0. So 
we can rewrite Eq. (|4.2p as 



At next-to-lcading order, we have 



sf 



a2g2C*(0)y 



(4.5) 



(4.6) 



where 



(0) 



K 



Inserting the leading order equation ()4.5p into (|4.6p . it gives 
This is the basic equation at next-to-lcading order. 



(4.7) 



(4.8) 
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B. Solution 



We can easily solve the leading order and next-to-leading order scalar field equations. The solution of Eq. (|4.5p is 
given by 



(0) 



MN) = Gf + mi{N -No) 



and the solution of Eq. (|4.8p is obtained as 







' I'M' 






/ dN"e 




J No ' 


J No 



sf 



Q2g2C*(0)y 



(4.9) 



(4.10) 



where A'o is an initial time and Cf and Dj represent the initial values of the scalar field and its time derivative. Note 
that the solution (°^(^/ satisfies the slow- roll conditions (j4.3p and (|4.4p if all the masses are small, 



Af2 = ^ < 1 . 



(4.11) 



Here let us show the time-independence of Sj, which is given by Eq. (|4.7p . From Eq. (|3.14p we have 
The leading order potential ^'^'>V is given by 



L I 



= Wexp 

where Cy is the initial value of the potential, 

Cy(x') = (°V(iVo) = Wexp 
Substituting the above solution into Eq. (|4.7p gives 



Cyexp M^{N~No) 



Y^miCf 



L I 



(4.12) 



(4.13) 



(4.14) 



. _ 3vW / e A 1 



D^{\ogCv) ~ -D'{\ogCv)D,{\ogCv) + D'{\ogCv)D'C^ 



mi 



R - Ud^C''' + 2I?*C"^ AC"^ ) + 2D'C'^jD,C' 



mi . 



(4.15) 



The time- independence of S'f is now clear since it is expressed solely in terms of C^, Cy, Cj and ^'^^^ij which arc all 
time-independent functions. Therefore we can rewrite the second order solution (j4.10p in a simpler manner as 



j2g2C'"(0)y 



where 



1 



Q2g2C'/'(0)y g2C* 



TV 



dN'a 



'„-3 



No 



N' 



dN" 



No 



(4.16) 



(4.17) 



Given the solution for the scalar field up to second order in gradient expansion, we can now obtain those for A', 
Aij and 7^ . At leading order, K is expressed in terms of the scalar field through Eq. (|4.12p as 



(0) 



K = V3(")y = V3CV cxp 



hl^{N-No) 



Because of the assumption (|2.2ip . ^^'^^ij and '■^^Aij arc trivial, 



(0) 



7y 



(4.18) 



(4.19) 
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At ncxt-to-leading order, Eq. p.l4|) becomes 



1 (2)y 1 



(4.20) 



and Eq. p.6p reduces to 



1 



a2e2C"^(o)if 



(4.21) 



where S:^ is independent of time. Specifically, 



TF 3 



-D, [D, (log Cv) j + -D, (log Cy)/?, (log Cv) 



3 r -1 T_F 



TF 



(4.22) 



Substituting the solution of the scalar field into Eq. (|4.20p . we obtain ^2)/^ as 



(2)K 1 



(o)K 6 



1 



5^, 



(4.23) 



where 



R - UD^C'f' + 2D'C'^'D^C* ) + D'CfD.Cf 



As for Aij, we obtain it by integrating Eq. (|4.2ip . 



(4.24) 



No 



Finally, '■^■'74^ is given by solving Eq. p. 71 



(4.25) 



N 



3{N'-No) 5-4 l-N 



6^ 



No 0. 



N' ^ 



(4.26) 



Before concluding this subsection, let us fix the remaining gauge degree of freedom on the uniform slicing 
mentioned in subsection III Dl We can make ("^AT spatially homogeneous on the initial slice by using this gauge degree 
of freedom. 



Cvix") = ("Vo = const. 



(4.27) 



where ''^•'Vq is a pure constant independent of both space and time. With this choice, we can substantially simplify 
the above expressions of the solution because the spatial derivative of Cy vanishes. 
To summarise, thus obtained solution is 



7ii 



Cf + mjiN-No) + ^Df 



K= V3(oyv\ l + -miD'f 
' 6 



2^^J.g3(N-Aro)j^(-^-) 2S',^J^(A^) 



C'' — 



a2e2C"*(o)T/ 



(4.28) 
(4.29) 

(4.30) 

(4.31) 
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where = (^Vo exp \^M'^{N — TVq)] and Cf, of, C^j and Cfj arc initial values of cj), dNcj), jij and Aij, respectively. 



The source functions Sf and are time-independent whose explicit forms are 



mi , 



Sf^ = 3 + D,C^D,C^ - D.DjC^ - -D,CfD,Cf 



TF 



and the fimctions J^, S , I a, and are given by 

r-W 1 /.AT' 

/Wo 



(o)-i/ ' 



1 1 (^)v 



1 r 



No 



(0)y 4 



/^(iV) = 0=^7(0)17 / dN' i 



JJN) = a^^°^V r dN' ^ I /"^ diV"^^ 



No 



(4.32) 
(4.33) 

(4.34) 



R - [iD'^C''' + 2AC"'' AC"^ j + ACf'ACf'J , (4.35) 

(4.36) 
(4.37) 
(4.38) 



C. Solution on the uniform K slice 

In this subsection, we derive the solution on the K gauge by applying a gauge transformation to the solution on 
the Af gauge. To do so, we first need to determine the generator of the gauge transformation between the two slices, 
N N = N + n{N,x'^) or conversely N + n{N,x^) = N. The nonlinear gauge transformation is discussed in detail 
in Appendix ICl 

As is clear from Eq. ()4.18p with the initial condition Cy = ^^■'Vb =const., Eq. ()4.27p . the uniformity of K on 
the uniform Af slices is kept in this model at leading order if it is chosen so on the initial slice. Thus no gauge 
transformation is necessary at leading order. 

At next-to-leading order, the gauge transformation of K is from Eq. (|C35P given by 



Since ^^''K vanishes on the uniform K slice, ^^-'n is determined from Eq. (|4.20p as 



2/1 



1 



Q2g2C'*(0)y 



The solution for tp is obtained from Eq. (IC33|) as 



2e 



3(N-N„) 



- 1 



1 



-s 



K 



where C"^ is given in p.2p and that for is same as 7^- on the M gauge from Eq. (|C34|1 . 



(4.39) 



(4.40) 



(4.41) 



(4.42) 



D. The curvature perturbation 

At next-to-leading order in gradient expansion, we have to extract the scalar component x from 7,^- to obtain the 
curvature perturbation [H. In linear theory, the variable x reduces to the traceless scalar-type component -ff^'", and 
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the curvature perturbation is given by T?.'"™ = {Hj^™ + _ff^'"/3)y, where ip —J- Hj^™ in the hnear hmit, and we have 
followed the notation of [i^l ■ 

Neglecting the contribution of gravitational waves and focusing on that of the scalar-type perturbations, we can 
assume that ^"■'7^ in the K gauge approaches the flat metric at late times when the adiabatic limit is reached, 



(0) 



7„- (5y (iV^O) 



(4.43) 



This condition completely fixes the remaining spatial gauge degrees of freedom, say, within the class of the time-slice- 
orthogonal threading. 

Under this condition, we rewrite 7^^ on the K gauge given in Eq. (j4.30p to manifest its time and spatial dependence, 



7., = -5,, + 2V3Cf^{x^) /^(O) + Cf^ix^) fHO) - 2V3Cf^{x^) /^(iV) - Gi^-(x') f^N) , 
where C7, = % from Eq. (|4.43p . and Cf, is given by 



(4.44) 



A 



TF 



The time-dependent functions / and are given by 



(4.45) 



(4.46) 



where the functions Ij and are defined in Eqs. ()4.36p and (|4.37p . respectively. 

Now we extract the scalar part from C/^ and C^j by using the definition of x given by Eq. ()3.20p . We note that 
this definition is unique in the sense that it reduces to the standard scalar part in the limit of linear theory and gives 
the 0{e^) correction unambiguously. 

The contribution of C/^ to x ni^Y be determined by evaluating Eq. p.4p on the initial slice with the help of 
Eqs. (IT^ and (|i3T|) . 



(4.47) 



where 



(4.48) 



Applying the definition of x in Eq. p.20p to the above, we find the contribution from is 



X^ = -y3(2)i^(iVo)(/^(0) - f\N)) + ^v^A-i9'(i?,^a,C/) (/^(O) - /^(7V)) 
As for the contribution of Cf^ to Xj we can make use of the relation. 



(4.49) 



-Ai? 



= d 



{■ 



(4.50) 



where the left-hand side is the Ricci scalar of e^*-^ 6ij . Thus the contribution from Cf^ is 



--R 



(/^(O) - ./'^(iV)) + ^A-ija^e-^'AaJ JgC^' (d.cfdjCfY^^ } (/^(O) - f^{N) 



(4.51) 



Adding both contributions together, we obtain 



x = x ■ 



(4.52) 



where x^ a-^d x^ are given, respectively, by Eqs. (|4.49p and (I4.5ip . 
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Finally summing all the contributions to the eurvature perturbation *H 
by Eqs. (|4.4ip and (|4.52[) . respectively, we obtain 



ipK + Xk/5, where ipK and xk are given 



1 



a2g2C'/.(0)y 



±(^)K{No)[f^{0)^f^{N) 



fHO) - f^{N)) + x!^A-'d^(^Dfd^Cf) (/^(O) - /^(iV) 



(4.53) 



What we need to know is the final value of U\ at sufficiently late times, N ^ (a ^ age We take A'o to be a 



time around which the scales relevant to cosmological observations crossed the Hubble horizon, hence iVo > 50. 
this case, at = 0, the curvature perturbation reduces to 



In 



D\k{N = 0) « ("^C^ + (2)C^ - ^_Df 



(4.54) 



The first term, ^^"IC^ represents the leading order curvature perturbation obtainable in the usual SN formalism, 
and the remaining terms represent 0{e^) contributions, the calculation of which is the main purpose of the beyond 
5N formalism. As a check, we have confirmed that the above result is consistent with the one obtained in linear 
perturbation theory with the same background solution. 

It should be noted that there is no contribution from xk at A'^ = by definition. However, this does not mean 
the evaluation of xk was meaningless. In general it can make an important contribution to 91 around the horizon 
crossing time N Nq, and by matching our 91 with that evolved from inside the horizon at ^ Nq, the final value 
9\{N = 0) is determined both by the values of tpK and xk and by their derivatives at iV ^ A^o- 

It may be noted that the last term proportional to Dj comes the time derivative of (f>i &t N = Nq, dN4>i{No) = 
of = 0{e^). Although this is generically small by construction, the contribution of the term itself can become large 
since it is proportional to mj/M'^ ^ 1/AI where Af^ ^ 1 as assumed in Eq. (|4.1ip . 



V. SUMMARY AND DISCUSSION 



In this paper, we developed a theory of nonlinear cosmological perturbations on superhorizon scales in the context of 
inflationary cosmology. We considered a multi-component scalar field with a general kinetic term and a general form 
of the potential. To discuss the superhorizon dynamics, we employed the ADM formalism and the spatial gradient 
expansion approach. 

Different from the single-field case, there is a difficulty in solving the equations in the multi-field case. At leading- 
order, the equations take the same form as those for the homogeneous and isotropic FLRW background with suitable 
identifications of variables. In particular, there are correspondences between the cosmic time in the background and 
the proper time along each comoving trajectory, t <^=> t, and the scale factor with the one defined locally ae^' a. 
This implies that given the background solution, the solution at leading-order in gradient expansion is automatically 
obtained. 

In the single-field case, one can show the coincidence among the comoving. uniform expansion, and synchronous 
slicings at leading order. This allows us to set the lapse function to unity, and replace the proper time by the cosmic 
time in the solution for the scalar field. Then the metric is expressed in terms of the scalar field easily, and the next- 
to-leading order equations can be solved straightforwardly because the space-time dependence of the source terms 
consisting of the leading order quantities is explicitly known. 

In cosmological perturbation theory, the most important quantity to be evaluated is the curvature perturbation 
on the comoving slices which is conserved on superhorizon scales after the universe has reached the adiabatic limit. 
This quantity accurate to ncxt-to-lcading order may be relatively easily obtained in the single-field case because of 
the above mentioned coincidence among the comoving, uniform expansion and synchronous slicings. 

On the other hand, in the multi-field case, such a coincidence between different slicings docs not hold. This implies 
the following. One can express the lapse function as a function of the scalar field, but the scalar field is also a function 
of the proper time. Thus one has the equation, 



/ t, 0(r) 



/ 



t, 



(5.1) 



To go beyond the leading order, one needs to solve this equation for a, but it seems almost impossible. 



16 



In this paper, we developed a formalism to go beyond the leading order which avoids the above problem. Namely, 
we first solve the field equations in a slicing in which the lapse function is trivial. The synchronous slicing is one of 
such slicings, but we adopt the uniform e-folding number slicing in which the time slices are chosen in such a way 
that the number of e-folds along each orbit orthogonal to the time slices, Af, is spatially homogeneous on each time 
slice. In other words, the uniform J\f slicing is a synchronous slicing if J\f is used as the time coordinate. 

In this slicing we can solve the equations to next-to-leading order without encountering the above mentioned 
problem. After the solution to next-to-leading order is obtained, we transform it to the one in the uniform expansion 
slicing which is known to be identical to the comoving slicing on superhorizon scales in the adiabatic limit. Thus the 
gauge transformation laws play an essential role in our formalism. We derived them which are accurate to next-to- 
leading order. Note that they are fully nonlinear in nature in the language of the standard perturbation approach. 
They are summarised in Appendix [Cj 

As a demonstration of our formalism, we considered an analytically solvable model and constructed the explicit 
form of the solution. Namely, we considered a multi-component canonical scalar field with exponential potential, 
= W eyip[J2jmj(j)j]. Following the general procedure discussed above, we first solved for the scalar field and 
the metric in the uniform Af slicing. Then using a remaining gauge degree of freedom of this slicing, we set the initial 
condition such that it coincides with the uniform expansion slicing at leading order. In this slicing with this initial 
condition, the next-to-leading order solution was straightforwardly found. Then we applied the derived nonlinear 
gauge transformation to it to obtain the solution in the uniform expansion slicing. Finally, from thus obtained spatial 
metric which takes the form, e'^^'ytj where jij is a unimodular metric, we constructed the generalised, nonlinear 
curvature perturbation defined by 



By inspecting the form of the final value of the curvature perturbation, we argued that the decaying modes of the 
scalar field may give rise to non-negligible contribution to the final, conserved curvature perturbation. To quantify 
the effect, it is necessary to match our solution on super-horizon scales with the one solved from sub-horizon scales. 
The evaluation of these next-to-leading order corrections to the spectrum as well as to the bispectrum of the curvature 
perturbation is left for future study. 
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In the single-field case, it can be shown that the four slicings, the comoving, uniform K, synchronous and uniform 
Af slicings, coincide with each other at leading order in gradient expansion. Here we demonstrate it. For simplicity, 
we consider a canonical scalar field, but the generalisation is straightforward. 

At leading order in gradient expansion, the Hamiltonian and momentum constraints are written as 




(5.2) 



Acknowledgments 



Appendix A: Coincidence of four slices in single scalar case 




(Al) 



(A2) 



where r is defined in terms of the lapse function as 




x'^— const. 




(A3) 
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Here let us choose the uniform K slicing, 

K{t,x^) ^ K{t) . (A4) 

Then Eq. (|A2[) immediately implies that the scalar field is uniform on this slice, 

-> = 0W- (A5) 

This shows the uniform K slice coincides with the comoving slice. 

Now since the left-hand side of Eq. (|A1[) as well as the potential term on the right-hand side is homogeneous, it 
follows that the kinetic term is homogeneous. This means 

drip = uniform — > a = a{t) . (-A-6) 

Therefore the uniform K slicing coincides with the synchronous slicing. Then the spatial homogeneity of both K and 
a implies that of the number of e- folds N by definition, Eq. (|2.36p . Thus we have shown the coincidence of these four 
slicings at leading order in gradient expansion. 

As we noticed in subsection III Dl there is a remaining gauge degree of freedom in the synchronous and uniform N 
slicings. In the above proof, this freedom is tacitly fixed since we took the uniform K slicing as a starting point. In 
other words, the uniform K (or comoving) slicing implies it is synchronous and uniform M . However, if one starts 
from the synchronous or uniform M slicing, we have to set the initial data such that K ov (p is uniform on the initial 
slice in order to show the coincidence. In fact, if we use this remaining degree of freedom to set a different condition 
on the initial data, for example '0 ~ 0, the synchronous or uniform M slicing will not coincide with the uniform K or 
comoving slicing. 

In passing let us also comment on the conservation of the curvature perturbation. From the definition of K, 
Eq. (|2.5p . one has 

dti}{t, x') = -H{t) + \a{t, x')K{t, x') . (A7) 

o 

This shows how the inhomogeneous part on the right-hand side contributes to the evolution of the curvature pertur- 
bation. As we have discussed above, there is no such inhomogeneities in the uniform K or comoving slicing in the 
single-scalar case. This is a proof of the conservation of the nonlinear curvature perturbation on super-horizon scales. 

On the other hand, in the case of multi-field, the corresponding scalar component on the right-hand side of the 
momentum constraint, Eq. (|A2[) . becomes uniform on the uniform K slice. However this scalar component docs not 
appear in Eq. (jAip in general, and hence one cannot show the homogeneity of the lapse function as in the single-field 
case. 

We note that this does not mean that the Hamiltonian and momentum constraints arc not related. Actually the 
leading order momentum constraint reduces to the spatial derivative of the Hamiltonian constraint in the slow-roll 
limit as shown in Appendix IE 2 al For a non-slow roll case, this relation between the Hamiltonian and momentum 
constraints holds only if we take into account the next-to- leading order corrections properly (sec Appendix IE 2 b[) . 



Appendix B: Uniform K slicing 



In this Appendix, we derive the basic equations on the uniform K slicing, K{t, x*) = 3H{t), and discuss a difhculty 
in constructing the solution at next-to-leading order in gradient expansion. 

In the uniform K slicing and time-slice-orthogonal threading, the constraint equations are given by 



7,2 p2V' 



6H^ = 2E , 



The evolution equation for K and Aij are 



R - AD^iP + 2DVAV' 



1 



drA, 



-iHAii - 



TF 



TF 



(Bl) 
(B2) 



(B3) 



(B4) 
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where 

and r is defined in terms of tlie iapse function as 

a{t' , x^) dt' , 



and the equation for the scalar field is 



,(ae^'P(ij)D'(t>' 



First, let us consider the leading order. The Hamiltonian constraint (|Bip reduces to 



(B5) 



(B6) 



(B7) 



(B8) 



hence ^'^^E is homogeneous on this slice. This means the solution for exactly coincides with the background 
solution Ebg- As for P, the leading order solution is easily constructed in terms of the background solution once it 
is known as discussed in Sec. Ill CI To be specific, if the background pressure is expressed as a function of t, 



PBG{t)^PBG t,Po{to) 



the leading order solution for P is given by 



The lapse function can be obtained from Eq. (jB3p as 



(0), 



dH{t) 



1 



dt (o)£:(t) + (o)P(t,a;*) 



(B9) 



(BIO) 



(Bll) 



Thus the lapse function becomes inhomogeneous if ^°'P is so. From Eq. (|2.5p . the evolution of -0 is related with the 
inhomogeneity of '■'^'a in this slicing as 



dt^ = H{a - 1) , 



(B12) 



which reflects the well-known fact that the curvature perturbation evolves if a non-adiabatic pressure exists. However 
there is a profound problem in Eq. (jBlip . 

Although it looks like Eq. (|Blip gives the solution for ("•'a, it is not so because the right-hand side depends also on 
'^'^^a. In reality, Eq. (jBlOp gives '■°^P as a function of t, which is given by the time integration of the lapse function. 
Therefore Eq (jBlip can be schematically expressed as 



a = f t,P{T) = f 



t,P 



idt 



(B13) 



In general it is very difficult to solve this equation, at least analytically. 

As clear from Eq. (|B3p . one needs to evaluate the spatial derivatives of a for example at next-to-leading order. But 
this will be almost impossible. This is the main problem in constructing a solution at next-to-leading order in gradient 
expansion on the uniform K slicing. This kind of problem arises also on other slicings except for the synchronous or 
uniform Af slicings. This is why the uniform J\f slicing is used in our formalism. 

As a closing remark of this Appendix, we emphasise that the above difficulty becomes actually problematic only at 
next-to-leading order. It is not a problem at leading order. This is because the solution of the "curvature perturbation" 
^^^ip is given by a time integration which depends only on the initial and final values of Af. From Eq. (|B12p we have 



x') - V(io, x') = / dt' - l)H = AA(t, x') - M{to,x') - N{t) + N{to) . 



(B14) 



Thus, at leading order in gradient expansion, as far as we are interested in the curvature perturbation there is no need 
to know the solution for '-'^'a. This is why the SN formalism works well despite the presence of the above problem. 
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Appendix C: Nonlinear gauge transformation 

We derive the gauge transformation rules for the metric, its derivative {K and Aij) and the scalar field. We consider 
a nonlinear gauge transformation from a coordinate system with vanishing shift vector = 0, to another coordinate 
system in which the new shift vector also vanishes, ^9* = 0. We note that once the time slicing is changed, the shift 
vector appears in the new slicing in general. So the spatial coordinates also need to be changed to eliminate thus 
appeared shift vector. 

We use the background e-folding number N as the time coordinate and define the temporal and spatial shift, n and 
L\ respectively, 

N ^ N = N + n{N, x'), 
x' ^ x' =x' +V{N,x'), 

or conversely 

N + n{N,S:') ^ N , (CI) 
x' +D{N,S:') = x\ (C2) 

Under the change of the coordinates, the line element should remain invariant, 

2 ^2 

= - rrZ^n ^N^ + a\N)e^'^j,jdx'dx^ = ^diV^ + a^{N)e^^^,jdx'dx^ . (C3) 

Equating the coefficients of diV^, dNdx^ , and dx^dx^ on both sides of the above, we obtain 

~9 9 

~\ 2 



(1 + %n) - {N)e^^^,,dfUd^U , (C4) 



= + df,n)(hn a^{N)e^'^^,,dfU{5\ + diL^) , (C5) 



1. Leading order gauge transformation 

First we derive the leading order gauge transformation. To begin with, we note the spatial shift U is 0{e) from 
Eq. (jCSp . since it is proportional to the spatial derivative of n, 

d^V ^ - 0{e) . (C7) 

Hence neglecting df^L'^ in Eq. (|C4p . we have 



This gives the leading order gauge transformation of the lapse function. 



(1 + 9^") a{N + n, x') + 0(e2) . (C9) 



H{N + n) 



The gauge transformation rules for the spatial metric components, ^ and 7y , are derived from Eq. (|C6[) . Taking 
the determinant of Eq. (jC6P and neglecting d^L^ and dih, we obtain 

a^e^V' = a^(N)e''^ + 0[^) . (CIO) 



20 



This gives the leading order gauge transformation of ip as 

i>{N, x') = ij{N, x')-n + 0(e2) 

= ij{N + n,S:')-r~i + 0(e^). (Cll) 

As is clear from Eq. (|C6[) . jij is invariant at this order, 

= 7^J(^ + n,^■') + 0(e2). (C12) 
Under the gauge transformation, the scalar field is essentially invariant except for the change of the arguments, 

4>{N,S:')^4){N,x'). (C13) 

Then, the leading order gauge transformation is 

(^(iV, i') ^ (P{N + i') + 0{e^) . (C14) 

Here we note that what we call a gauge transformation should perhaps be called a coordinate transformation as far 
as the change of time slicing is concerned. For example, for a scalar quantity Q, the gauge transformation under the 
temporal shift (t — > t + T) is given by 

5Q(t, x') = 5Q{t, x') - TdtQo , (C15) 

which is derived from the invariance of Q as a scalar under the time coordinate transformation, 

Q{lx')^Q{t,x'). (C16) 

In the context of gradient expansion, we define the nonlinear gauge transformation of a scalar quantity by Eq. (jC16P 
not by Eq. (|C15[) . since the gradient expansion is a completely nonlinear approach. If we were to expand perturbatively 
the left-hand side of Eq. (|C16[) . we would obtain an infinite number of terms, which is definitely something to be 
avoided. On the other hand, as we see shortly below, to the accuracy of our interest, i.e. to 0{e^), the spatial 
coordinate transformation can be linearised. Hence we interpret it as the standard gauge transformation and compare 
the quantities at the same coordinate values. 

Once we have the gauge transformation rules for the metric components, we can readily obtain those for the extrinsic 
curvature components, K and Aij. Since Aij vanishes at leading order, we only have to consider the transformation 
of the expansion K. From its definition, Eq. (|2.5|) . and using the transformation rules for the metric components 
derived above, we find it is invariant at leading order, 

k{N, = K{N + ("^n, x') + 0{e^) . (C17) 

Again, as discussed around Eq. (jC16[) . we note that K is invariant only in the sense of the nonlinear gauge transfor- 
mation we have defined. 



2. Next-to-leading order gauge transformation 

Now we derive the gauge transformation at next-to-leading order. To do so, we first have to fix the spatial coordinate 
transformation. It is determined by Eq. (jCSp . which results from the requirement that the shift vector should vanish 
to O(e^). In terms of fi, the gauge transformation generator L* is given as 



a?{N)e-^^m{N) 
' dN — 



where P is a time-independent spatial vector which represents the remaining gauge degrees of freedom in the spatial 
coordinates in the time-slicc-orthogonal threading. 
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To consider the gauge transformation at next-to-leading order, wc expand the nonlinear gauge transformation 
generator fi in Eq. (jCip to the leading order term, the next-to-leading order term and so on, 



(C19) 



As for the spatial shift L*, it is unnecessary to expand it, since it is 0{e) already and the ncxt-to- leading order term 
in U is 0{€^) which only affects the gauge transformation at 0{e^). 

First we consider the metric components. From Eq. (|C4p . the next-to- leading order gauge transformation for the 
lapse function is determined as 

_ LX^«)4^(l + a^(°)n) + Oi.^) , (C20) 

2a2((0)7v)e2'Ai73((0)^) ' 'J 'h{(o)N)^ ^ ' ^ ' ^ ' 

where di — d/dx^ and '^°'iV(A^, a;*) = + '"'ri(A, x*). Note that all the quantities on the right-hand side are those 
evaluated at N ~ N and a;* = 5:*. Taking the determinant of Eq. (|C6p . 



we obtain the gauge transformation for ^^^V'l 



^ ^ 3 ' 6a2((0)Ar)e2V'H2((0)^) 



Then using Eq. (jC6p again, the gauge transformation for is obtained as 



a2((o)^)e2<AiJ-2((o)7v) 



Next we consider the extrinsic curvature. At next-to-leading order, Eq. (|2.5p gives 

(1 + Ofr^^hija \ a 1 



^(0)^5 (0)^ 



2(l + a^(0)n)a2((0)A)e2^i7((0)A) ' 2aa2((o) A)iJ((")A) Ve^''' 



3g(WA) 
(l + 9^(o)n)a 

iJ((o)A)a7*J' 



2(l + 9^(0)n)e2'A ^ \^a2((o) A)i72((o)7V) 
And also Eq. (|2.16p determines the gauge transformation for Aij as 

i?((0)A) 



A,,(A,i') = A 



2(l + aj^(o)n)a 



2(l + aj^(°)n)a 



TV 



a2((0)iV)e2'Ai/2((0)7v) 



Finally, the next-to-leading order gauge transformation for the scalar field is given from Eq. (jC13p as 



(C21) 



,^yja^(o)^5 (o)~ _,_0(g4)^ (C22) 



(^d'^°^ndf'>h - i7"a^,(°)na/o)n7„-^ + 0(e4) . (C23) 



(C24) 



(C25) 



(C26) 
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To conclude this Appendix, let us summarise the nonlinear gauge transformation derived above. The leading order 
transformation rules are 



a{N, = ^/^j (1 + % Wn) a + 0{e^) , 
7,j(iV,i*)=7,,+0(e2), 

The next-to-leading order transformation rules are 



H{(.o)N) 



2a2((0)Ar)e2>/.iy3((0)^) ' » J ^ * ^ 77((o)Ar) ^ ^ ^ ^ ^ 



6a2((o)7v)e2'Aif2((o)^) 



7.,(iV, i') = + (2)na^7,^. + i'a^7., + lok&iL" + 7,^.9^^' - ^d^L^, 



a2((o)iV)e2V'i/2((o)7v) 



and 



(2) 



(1 + df,(^'>n)a \ 



N " 



(2) a a^(2)n 



a(0)^n (0)^ 



2{1 + dff(°)n)a^{(")N)e^^H{(0)N) ^ 2aa2((o)7v)iJ((o)7V) ^ U^"'^ 



3g(WjV) 
(l + %(0)n)a 

i/((")iV)a7*J' 



^(o)^a (0)^ \ 



■(l + c>^(")n)e2'A ^ ^a2((o)^)7j2((o)^)y ^ ^ > 



4- (iV, 5^) = A,, - + jj.djdj^L'^ + i.^kd-dfL'^) 



TF 



2(l + 9^(")n)a^ 



.a2((0)7V)e2'Ai72((0)7v) 



TF 



where is given by 



No 

((0)iV)e2'Ai/2((0)^) 



(C27) 

(C28) 
(C29) 
(C30) 
(C31) 



(C32) 
(C33) 

(C34) 



(C35) 



(C36) 
(C37) 

(C38) 



Appendix D: Consistency checking 

In this Appendix, we apply our formalism to the case of a single scalar field and verify the consistency. One way to 
do this is just to compare the obtained solutions in this formalism with the results obtained previously, for example, 
in j28j . Here, however, we take another, direct way. Namely, we consider the solution in the comoving gauge and 
transform it to the one in the Af gauge. Then we insert thus derived solutions for the metric and the scalar field into 
the basic equations in the M gauge and show that they are satisfied. 
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1. Solution in the uniform M gauge 

First we construct the solution in the M gauge by transforming it from the one in the comoving gauge. As we have 
seen in Appendix \^ the leading order solution in the comoving gauge is 

= /„(Ar), (o)^e=C(x'), ^''h^Jc^C,J{x'), K, = fK{N), 0, = /^(iV). (Dl) 

Here fa, fx and are functions of only time and have no spatial dependence. On the other hand, C and Cij are 
functions of only spatial coordinates and have no time dependence. Because of the definition of K, Eq. (|2.5p . the 
functions fa are /„ are related as 

fK = ^ . (D2) 

Ja 



At next-to-leading order, however, all the quantities listed in Eq. (|D1[) become space-time dependent except for the 
scalar field ^^^<^, which is by definition a function of only time. Hence without loss of generality, we may absorb it in 
the leading order scalar field and set '■^V = 0. 

To perform the transformation from the comoving slicing to the uniform M slicing, N = N + h{N, i*), one needs to 
determine n. In the uniform Af slicing, tp is time-independent by definition, Eq. (|2.36p . Considering the transformation 
of tjj, this gives at leading order, 

^A^(F) -C(a;')-(0)n + O(e2) ^ = C(F) - + ©(e^) . (D3) 

At the next-leading-order Eq. (|C33[) gives 

' 3 ' 6a2((0)Ar)e2CiJ2((0)^) ^ > 

(2)n = (2)^^ + L'd-fi + -chV - - , ' ^ + 0(e^) , (D4) 

^ ^ 3 ' 6a2((o)Ar)e2Cij2((o)7v) ^ ^' ^ ^ 



where is given by Eq. (|C38p , which is derived from the time-slice-orthogonal threading condition, and the fact that 
'•^■'V'A/' = is used since it may be absorbed into '•^-'i/'AA- Note that from Eq. (jD3[) . '■°^ri is independent of time which 
significantly simplifies the analysis below. 

Performing the gauge transformation, we we obtain at leading order, 

«AA(iV, i') = -^i^/^WiV) + 0{e') , (D5) 

7,,Ar(7V,i*)=C„+0(e2), (D6) 

KMN,i')=fKi^"^N) + 0{e'), (D7) 

MN,i')=f4.{^°^N)+0{e^), (D8) 

and at the next-to-leading order, 

' H{^°)N)\ U U Hi^N) J 2 a2((0)^)e2Ci/3((0)7v) ^ ^' 

(D9) 



J a 



a2((o)Af)e2Cff2((o)7v) 



(DIO) 
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and 

(2) 



2e 



2C 



2/a 



a2((o)iV)e2CiJ((o)iv) 



a2((o)Ar)if2((o)7v) 



TF 



J Oi 



a2((o)Ar)H2((o)7v) 



where di — d jdx^ and U is given explicitly by 



o2C 



dN'- 



J a 



I No ^ 

After some manipulations, Eqs. (|D9|) . (|D11|) and (|D12|) may be re-expressed as 



^ ~ ? - d/'^n - ('^nd^ log /^((°) TV) - ^ 



9 D^^°^fiD^S°'^n 



(0) f 



2a2((o)Ar)e2C/2.((o)7v) 



(Dll) 

+ 0(e4), (D12) 

(D13) 

©(e^) . (D14) 

(D15) 



- - (2)~5^^^((0)^) + 3 ( _ + log/,^((")7V) 



n2 (0)^ I r)j(o)^ n 
3 , i , (D16) 



a2((o)7V)e2C/^((o)7v) a2((o)iV)e2C/^((o)7V) 



-AijJ\f Ai 



a2((0)iV)e2C/;^((0)iV) 



j((o)7v)e2C/^((o)Ar) 



^(l-dj,\ogfKrN))[D,^'^nDf^^ 



TF 



TF 



(D17) 



where is a covariant derivative with respect to ^^'^jijj\f (or equivalently to Ctj). These expressions are convenient 
for later use. 



2. Consistency check for Aij 



Here we explicitly show that the solution for Aij obtained in the previous subsection satisfies the evolution equation 
in the Af gauge. This verifies the consistency of the derived gauge transformation. 
In the J\f gauge, the evolution equations for Aij is 



d,rA 



3A 



3 



TF 



n TF 



while in the comoving gauge it is 

dNAijc = 3A 



■ [r,j + DiCDf - D.DjC 



TF 



'''^ a^N)e^CfKiN) 
We rewrite the second line inside the brackets of Eq. (jPlSP as 

(2nd line) = (% log /k) {^°^N) [D.D/'>n - (o.'^^'^nD.^u + Df^nD-i^^^ 

d%log fK - (9^1og/A')' - 2%log/^,l (W7V)i?jWni?.(0)n, 



(D18) 



(D19) 



(D20) 
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where we have used the leading order Einstein equation, 

By subtracting Eq. (|DT9|) from Eq. ([DTsI) . we find with the aide of Eq. ((D20)) . 



(D21) 



lA/" 



^ g_f + %log/K((°)^) 



a2((o)Ar)e2C/^((o)Ar) 
3 



TF 



„2((o)^)e2Cj^((o)7v) 



(D22) 



where Q I ^ = QaT - Qc- 

On the other hand, taking the time derivative of Eq. (|D17|) . which we have obtained by the gauge transformation, 
we obtain 



3A, 



a2((o)Ar)e2C/^((o)Ar) 
3 



i 1 



TF 



a2((o)7v)e2C/^((o)7v) 
Comparing Eqs. ((D22|) and ((D23)) . we see the precise coincidence between the two. 



TF 



(D23) 



3. Consistency check for K 



In the 7V^ gauge, the evolution equation for K is 

iH{N),^ , ,2 3 



2aM 



■(9jv0A^) - 



2a2(iV)e2'/'AAAV 



and that in the comoving gauge is 



2ac ' 2a^N)e^^^fKiN) 

Again we rewrite the second line inside the brackets of Eq. (jD24[) as 



R - AD^C + 2D'CD.,C 



(D24) 



(D25) 



(2nd line) ^ \^{dj;; log JkY ~ dllog Jk] N)D'^^>^hD-'^"^h - dj^ log fxi^^^ , (D26) 

where we have used Eq. (|D2T|) . 

By subtracting Eq. ((D25|) from Eq. ((D24)) . we find with the aide of Eqs. ((DTs)) and ((D26)) . 



%log/A-((°)iV)-2 
'a2((0)7V)e2C/^((0)7V) 



a2((o)iV)e2C/^((o)iV) 



1 + :^aA.log/K((")7V) + 4log/,,(WiV) - (%log/K)'((°'iV) 



(D27) 



On the other hand, taking the time derivative of Eq. (|D16p , which has been obtained by the gauge transformation, 
we obtain 



(^^^ndMfKiN)) +3^ 



a2((o)7V)e2C/^((o)iV) 



dj^\ogfK{^^)N)-2 
2((0)iV)e2C/^((0)7V) 



1 + -a^iog/K((")7v) + 4iog/K((")iv) - (a^iog/K)'((")iv) 



(D28) 



We see that Eqs. (jD27p and (jD28|) coincide precisely with each other. 



26 



Appendix E: Note on constraints 



In general relativity, one cannot freely choose the initial values of the metric and matter field variables. They must 
be chosen so as to satisfy the Hamiltonian and momentum constraints. In this Appendix we make some comments 
on these constraint equations in the context of the spatial gradient gradient expansion. 

Let us consider M scalar fields minimally coupled with Einstein gravity and count the number of degrees of freedom 
in this system. In the ADM (Hamiltonian) formalism, the dynamical geometrical degrees of freedom are six spatial 
metric components {tp, jij) and their time derivatives {K^ ^ij)- In addition, we have 2Ai dynamical degrees of freedom 
associated with the scalar fields and their derivatives. To summarise, 

C^:l, C^^:l, C,^ : 5 , Cf : M , of : M , (El) 

that is, we have 12 + 2Ai (= 6 + 6 + A^+A^) degrees of freedom. The Hamiltonian and momentum constraint 
equations give four constraints among them. And there are also four gauge degrees of freedom. Thus there are 
4 + 2Ai (= 12 + 2A4 — (4 + 4)) independent physical degrees of freedom left, which represent 2x2 gravitational wave 
degrees of freedom and 7W x 2 scalar field degrees of freedom. 

Let us explicitly check how the above counting works in the Af gauge, where we take the uniform Af slicing and 
the time-slice-orthogonal threading. The 2M scalar field degrees are represented by Cf and Df. The gravitational 
wave degrees of freedom are contained in among 5 degrees of freedom in C^^ and C:^, respectively, 2 of each describe 
the gravitational degrees of freedom. The remaining 3 degrees of freedom in C^j are fixed by the threading condition 

= together with purely spatial coordinate transformation degrees of freedom (see e.g., Eq. (|4.43|) in the example of 
a canonical scalar field), and those 3 in are fixed by the momentum constraint (j4.47|) . The Hamiltonian constraint 
is used to determine . 

Now we are left with C^. This may be regarded as the remaining gauge degree of freedom in the Af slicing, since 
the slicing condition implies dttp = 0. But once the scalar field configuration on the initial slice is fixed, it should 
not be freely specifiable, since if it were it would contradict with the total number of the physical degrees of freedom 
counted above for the general case. In subsection IE 11 we find there is indeed a constraint equation that must be 
satisfied by C^' . Thus we recover the total number of the true physical degrees of freedom correctly. 

Recently the validity of the 6N formalism was studied by Sugiyama et al. [s^l, and they found the violation of 
the momentum constraint in a non-slow-roll inflation case. In subsection IE 21 we show how the consistency with the 
momentum constraint is recovered in the gradient expansion. 



1. "Hidden" Hamiltonian constraint 



In general the Hamiltonian constraint gives a non-trivial relation among the initial values of the metric components, 
their time derivatives and the matter field variables. However, in Sec. IIVI we have used the Hamiltonian constraint 
to solve for K. Thus it seems it will not give any more constraint. Then where is a constraint corresponding to the 
Hamiltonian constraint? The answer is that it is hided in the evolution equation for K. 

In order to see it explicitly, let us consider the canonical scalar field model discussed in Sec. IIVI We first solve the 
evolution equation for K, Eq. p.5|) . and then compare its solution with the solution obtained form the Hamiltonian 
constraint. 

Let us consider the evolution equation for '^•'/v in the J\f gauge. Using the leading order solution and given by 
Eq. (|4l0|) . we find 



dN' 



io)v 



1 



2a2e2C"*(o)y 



R - {aD^C"^' + 2D*C"'' AC"^) - D'CfD.Cf 



(E2) 



This is easily integrated to give 



V3(o)y 



[N-No) 



miSp^jN) 



R - iD^Cf' + 2D*C"^ AC"^ - D^CfDiCt / r i 



2a2e2CV'(o)v^ 



2(0)y / ' 



(E3) 
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where the suffix tmp indicates that it is obtained from the temporal evolution equation. 

Now let us carefully examine the solution (|E3|) in comparison with that obtained in the text, Eq. (|4.29p . First we 
compare the initial values. This determines the initial value of '■^•'-K^tmp as 



R - Uo^C'f' + 2D'C"f'DiC'f') + D'CfD,Cf 



Next we subtract the solution (|E3p together with thus obtained initial data (jE4p from the solution (|429)) . 



(E4) 



tmp 



niiSf 



1 f 1 



1 



2a2g2C'^'(0)y 

After integration parts and some manipulations, we obtain 



(2)A--(2)A- 



tmp 



1 



2a2e2C'*(0)y 1 eC* 



D'CfD.Cf 



R - [AD^C^' + 2L>'C"^£'iC"'' ) + D'CfD,Cf 

1 



(E5) 



mi 



D 



AP + 6 
6 



Af2 



R 



(E6) 



The right-hand side of the above equation must vanish at all times. This demands that the term inside the square 
brackets should vanish. 



6 



R - ( AD^C'f' + 2D'C'''DiC'^' 



= 0. 



(E7) 



Apparently this gives a non-trivial constraint among initial data, and Cf. This is the "hidden" constraint 
corresponding to the Hamiltonian constraint. 

2. Consistency with the momentum constraint 

Below we show that the momentum constraint is automatically satisfied if the Hamiltonian constraint as well as 
the scalar field equation is satisfied. In the case of slow-roll infiation, the leading order Hamiltonian constraint is 
sufficient to show it. For general (non-slow-roll) inflation, one needs the next-leading order Hamiltonian constraint. 

a. Slow-roll case 



Under the slow-roll approximation, the Hamiltonian and momentum constraints reduce respectively to 

2 

— J 

3 



2 2K , . 

-d,K = —P(ij)dN(l>'d,<^\ 



(E8) 
(E9) 



and the scalar held equation becomes 

Now we multiply both sides of the momentum constraint (jE9p by 2A', and substitute the scalar held equation (jE10[) 
into the right-hand side of it. This gives 



(ElO) 



(left - handside) = '^KdiK = d, Q^^^ . 
(right - handside) = -2PA4>^ = d,{-2P) . 



(Ell) 
(E12) 
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One can easily see that these are merely a spatial derivative of the Hamiltonian constraint (jESI) . Thus we conclude 
that, at leading order in gradient expansion, the momentum constraint is automatically satisfied if the Hamiltonian 
constraint and the scalar field equations are satisfied. 



Non slow-roll case 



Here wc show that in general the leading order momentum constraint is automatically satisfied once the Hamiltonian 
constraint and the energy conservation equations are satisfied to next-leading order in gradient expansion. 
From Eqs. p.f 3p and (|2.14|) with Aij = ©(e^), the Hamiltonian and momentum constraint equations are 



1 



R 



The energy conservation law, n'^'V,^T'^'^ = 0, in the Af gauge is 



OnE = 3(E + P) 



-diK = J, 



2E + 0{e'^) 



(E13) 
(E14) 

(E15) 



Now let us multiply the evolution equation for K, Eq. p.5|) . by 2K/3 and subtract Eq. (jEf Sp from it. This gives 



On ( ^K' 



E 



I 



R- (4D^i; + 2D'i;D,iP 



2K 



K 



K 



3K 

Q2g3V 



K- 



Substituting the Hamiltonian constraint (|E13[) into the left-hand side of the above equation, we find 



OnIIk'-E]^ -d. 



Comparing the above two equations, wc obtain 



R- Ud^v + s^VA^ 







R- (4D^ip + 2D'ipD,i; 



(E16) 



(E17) 



-diK - J, 



(E18) 



Again one can see that the leading order momentum constraint (jE14p holds automatically if the next-to-leading order 
Hamiltonian constraint (|E13|) holds. 

What docs this mean? In the gradient expansion, wc have assumed Aij does not contribute to the leading order 
dynamics. This corresponds to neglecting the adiabatic decaying mode in linear theory. In a single field model, it is 
well known this decaying mode appears only at the next-leading order both in linear theory j3l| and in non-linear 
theory [2^ [2^ . We also mention the work in which they discussed the behavior of decaying modes in different 
choices of gauge. The absence of the decaying mode at leading order in gradient expansion should hold also in the 
case of multi-field infiation, at least as long as the background homogeneous solution is stable against a homogeneous 
but anisotropic perturbation. 

To summarise, the point is that the initial data for the scalar field and its time derivatives arc not freely specifiable 
in general, and in particular we need to take into account the next-leading order terms in the Hamiltonian constraint 
for the non-slow-roll case. Once we take into account the next- leading order Hamiltonian constraint, the leading order 
momentum constraint is automatically satisfied. As for the ncxt-to-leading order momentum constraint, it constrains 
the initial value of Aij as in Eq. (|4.47p . 
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